On Solutions for Two Waves with 
Periodic Coupling 

By S. E. MILLER 

(Manuscript received January 25, 1968) 

An exact solution jor the coupling effects between two waves with a 
■particular complex periodic coupling junction is presented; the particular 
coupling junction gives the same wave interactions as constant coupling but 
at a translated value oj differential phase constant. A transjormation is 
given which permits known theory jor constant coupling to be applied to the 
periodic coupling case. 

Approximate solutions are given jor periodically reversed coupling 
(sinusoidal or square wave) between two waves, and calculations are pre- 
sented which indicate the solutions are valid jor arbitrarily long coupling 
regions or arbitrarily large integrated coupling strengths. The region oj 
validity jor earlier perturbation theory is defined and proved to include the 
cases oj interest jor multimode circular electric waveguides. 

I. INTRODUCTION 

This paper describes some solutions for two waves with periodic 
coupling. Coupled waves have been important in a wide variety of 
communication devices: transmission lines, directional couplers, am- 
plifiers, and in describing mode interchange phenomena generally. 1 
Multimode transmission lines have been advantageously described 
through coupled wave equations, and a particular situation of im- 
portance exists in the circular electric waveguide. 

As first shown by H. E. Rowe and W. D. Warters, 2 periodic 
straightness variations cause periodically reversed coupling from the 
circular electric wave to several other waves, and this interaction 
results in the most difficult tolerances on the fabrication and instal- 
lation of the waveguide itself. Publications by H. E. Rowe and 
W. D. Warters have provided a comprehensive understanding of the 
fundamentals involved and have given explicit expressions based on 
perturbation theory for calculating the loss versus frequency varia- 
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tion resulting from such periodic mode conversion. 3 A subsequent 
publication by D. T. Young* has indicated how the approximate 
solution of Rowe and Warters may be transformed mathematically 
to show explicitly the way differential attenuation smooths out the 
mode-coupling effects on the loss versus frequency characteristic. 
Young's solution depends upon a valid perturbation solution. A 
discussion of the accuracy of the perturbation solution is given by 
H. E. Rowe; 5 when the differential attenuation is too small the per- 
turbation solution breaks down and it is of interest to know exactly 
where and quantitatively how this occurs. 

We show in this paper that a certain periodic coupling function 
has exactly the same effect on waves of unequal phase constant as 
uniform coupling between waves of identical phase constant. A trans- 
formation is given to allow the use of earlier theory for periodically 
coupled waves. 

Also presented here is an approximate solution for the periodic 
coupling distribution sketched in Fig. 2 valid for any value of dif- 
ferential attenuation. It is true that known solutions for uniform 
coupling, as in Fig. 1, can be applied to Fig. 2 by simplying solving 
for the output values at x = A,„/2 and using these as the input bound- 
ary conditions for the transmission region starting at x = \ m /2. The 
resulting exact expression representing conditions at x = k m can be 
expressed as a matrix and raised to the n th power to represent the 
solution at x = nA™. We seek here a simpler form of expression in 
which the functional interrelations can be visualized without exten- 
sive numerical calculations. 

II. EXACT SOLUTION FOR TWO PERIODICALLY COUPLED WAVES 

We start with the following equations for two coupled waves : 

f^M = - 7 i#, +c ai (z)E i (1) 

j z E 2 (z) =c l2 {z)E i -y 2 E 2 (2) 

c(z) 



Fig. 1 — Constant coupling. 
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Fig. 2 — Sq»iare-wave coupling. 

in which yi and y 2 are the complex propagation constants and (fa 
and C21 are coupling functions. Appendix A shows that the coupling 
functions 

c ai = jc„ exp (jk e z) = jc p [cos k c z + j sin M] (3) 

and 

c I2 - jc p exp {-jk<z) = ic p [cos M - j sin M] (4) 

give the same solutions previously found 8 for k c = provided that 

k c = (0, - /SO = A/3 (5) 

where c p is a constant and 

7 m = a m + J/3 m . (6) 

When (5) holds, complete transfer of power between waves can 

occur. When (5) does not hold, the resulting wave interactions can 

be calculated using previously developed theory for k c = and 

substituting AB for A/3 in the fc c = solutions, where 

AB = A/3 + A-, (7) 

III. PHYSICAL REALIZATION OF IDEAL PERIODICALLY COUPLED WAVES 

We describe here a physical realization of waves coupled according 
to equations (3) and (4), and cite an advantage in mode selective direc- 
tional couplers. 

Figure 3 shows a mode-selective coupler between TE° of rectangular 
guide and TB& of round guide. The thin dielectric lining is used to break 
the TE^ — TMf t degeneracy. The longitudinal magnetic intensity h% 
of the TE& wave is coupled to the longitudinal magnetic intensity h° 
of TE° in the off-axis longitudinal slots and is also coupled to the trans- 
verse magnetic intensity h° of TE° in the 45° slots on-axis. In each case 
the magnitude of the coupling is set by the length and width of the slot. 
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Fig. 3 — TE°i -TE ° i coupler using unequal phase constants in the rectangular and 
round waveguides. 



The phase reversal of the longitudinal h° coupling is accomplished 
by reversing the slot position relative to the centerline, and the phase 
reversal of the transverse h° coupling is accomplished by reversing the 
slant angle of the slot. One set of slots represents the sine term. It is 
not necessary that any particular fraction of X m be used in a coupler, 
since the coupling is always in the same phase relative to the desired 
waves in the two guides. To accomplish the desired mode selectivity 
the A/3 between the TE° Q and TE& waves is made equal to k e = 2ir/\ m 
as in equation (5). 

An advantage of a coupler of this form, compared with one in which 
constant coupling is used with A/3 = 0, is that the waveguides can have 
the standard dimensions set by other considerations. 

Other illustrations of useful coupling between waves of unequal phase 
constants will be given in another paper which the author is preparing. 



IV. SQUARE-WAVE OR SINUSOIDAL COUPLING 

We present here the results of Appendices B and C which discuss 
approximate solutions for the cases in which the coupling is a square 
wave as in Fig. 2, or the corresponding sinusoidal 



c(z) = sin (2irz/\ m ). 



(8) 



Both solutions are expressed in the form and notation of a previous 
publication 6 giving the solution when the coupling is constant, and 
the boundary conditions Ei(0) =1.0 and E 2 (0) = are impressed: 



E y (z) = e- y, '{Ae T ' ! + Be'*'! 



EM - f^= W" - e"\ 



(9) 
(10) 
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(12) 

(13)* 
(14) 



Aa = a x — a 2 

A7* = Aa + i A/3. . 

The wave interactions are described by the above equations provided 
the following values of c* and A/3* given in Table I are used. 

Table I — Values of c* and A/3 



Coupling function 
(see Fig. 2) 


C 


A0. 


Constant = c 


C 

2 

- c 

7T 

c 
2 


A/3 = (/3i - ft) 


Square-wave of 
magnitude c and 
period X m 

C Sin (If) 


«-£>M*y 



These solutions, equations (9) and (10), have been obtained in Appen- 
dix B by relating the rate of transfer of power (that is, transfer over a 
short length interval) for the periodicaUy reversed coupling to that for 
constant coupling, and noting the effective value of coupling c* and 
effective differential phase constant A/3* . The solutions are correct 
for z equal to an integral multiple of X./2, and may be in error by less 
than approximately 0.2cX m /ir at intermediate values of z. 

There probably should be a correction factor in A/3* for sinusoidal 

*r, corresponds to the + sign and r 3 corresponds to the - sign. 
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coupling similar to the radical shown in Table I for the square wave 
coupling; the work done thus far has not defined what it should be, but 
for small cX m the radical is negligible for most purposes. 

The values of A/3* in Table I bear marked resemblance to equa- 
tion (7) ; however, because the simple sinusoid or square wave cou- 
pling phase is two-valued (versus z) instead of continually progressing 
(equations 3 and 4) to provide coupling continuously in step with 
the phase changes of waves 1 and 2, there are local maxima* in 
coupling effects at other values of A/3 for sinusoidal or square wave 
coupling. Appendix B shows that the wave interaction effects are 
properly described for square wave coupling in the regions near 

A/3X. = 2*p (15) 

p = 1, 3, 5, ••• 
by the transformations 



Kpir 



A/3, = A/3-pf A /l -(?£=) (16) 



c* = *. (17) 

* pw 

Because c, drops off rapidly with increasing p the corresponding wave 
interaction effects drop off also. 

4.1 Numerical Comparison of Approximate and Exact Solutions 

A few calculations have been made to find quantitatively the 
error resulting from the approximations made in equations (9) and 
(10) for coupling as in Fig. 2. An "exact" solution is obtained by 
using exact uniform coupling theory on each interval of 0.5 a™, the 
output of one interval being taken as the input to the next interval. 

We take first the simplest case, 

Aa = 0. 
Then equation (9) becomes 

• jcos [ V~ I cz] - i^jd sin-[ V" I cz]\ 1 s 



* Rowe and Warters noted this in their work recorded in Ref. 3. 
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(20) 



and A/3* is as given in Table I for square wave coupling. We notice in 
passing that for A/3* = 



(21) 
(22) 



Ei = cos \-cz 



Ei = i sin ( - cz 



and the power exchanges completely back and forth between the 
two waves as a function of z; this is of course identical to the be- 
havior in uniformly coupled waves, but with a modified period given 
by the 2/tt factor. 

We take for the first numerical comparison the condition 

2 7T 

-cz = 5 

so that at A/3* = 0, E x — and \E 2 \ = 1.0. The additional specific 
numbers used are 
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Fig. 4 — Undriven-wave amplitude versus distance for lowest order square- 
wave coupling. 
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X m = 2 feet 

z = 1000 feet 

c = 2.46740 X 10 -3 ft" 1 . 

Fig. 6 shows a plot of the loss 20 log | E x | (labelled | Aa/c | = 0) versus 
A/3 + X m /2 and Table II shows the comparison of the exact versus approxi- 
mate calculations, the latter obtained from equation (9) and associated 
expressions. It may be kept in mind that A/3 and A/3,,, are inversely pro- 
portional to frequency in many cases of interest, so that Fig. 6 is a loss 
versus frequency plot associated with a particular periodic coupling 
component. The radical in the expression for A/3.,, was ignored in the 
comparison of Table II and this might account for the consistent positive 
difference (approximate — exact) for A/3X m /2 greater than it. Aside from 
the pole at A/3X m /2 = t, the two calculations agree to better than 1 
percent in dB even when the loss is a few tenths of a dB. 

Figs. 6 and 7 and Tables III, IV, and V show similar comparisons 
for Aa = — | c |, A« = — 5 | c |, and Aa = —50 | c |. Excellent agree- 
ment is obtained in all cases. 

Table VI shows a comparison of the phase angle on E x , computed by 
the two approaches. For the points shown and for the other points 
(not shown) corresponding to the amplitude values of Tables II and V 
the agreement is excellent. Fig. 8 shows a plot of the phase, where odd 
symmetry about A/3 + = is understood. 

A check has also been made on the accuracy of equation (9) in the 
region near A/3X m = 6tt, corresponding to p = 3 in equations (15) through 
(17). The same parameters were used as in the calculations for Figs. 
6 and 7. The results are plotted in Figs. 9 and 10 which represent both 
the approximate calculation from equation (9) and the exact calcula- 
tion. The differences are on the same order as given in Table II, and 
are too small to show in the figures. Figures 9 and 10 may be compared 
directly with Figs. 6 and 7 to see the "third harmonic" loss (labelled 
E\ in Figs. 9 and 10) in relation to the "fundamental" loss, Figs. 6 
and 7. 

4.2 Interpretation and Further Simplification 

Consider first the shape of the loss versus A/3 + (equivalent to loss 
versus frequency) curves. In the limit 

! Aa 
— » 1 (23) 

c 
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Fig. 5 — Undriven-wave amplitude versus distance for "third harmonic" in 
square-wave coupling. 



Table II — Exact and Approximate Calculations 

OF Ei FOR | Aa/c I = 



A»fc\ 


20 log 


1 A| 


Approximate 


^Kl) 


Approximate 


Exact 


-exact 


w - 0.0157079 


- . 1662db 


- 0.1680db 


+ 0.0018db 


w - 0.01413716 


- 0.1389 


- 0.1403 


+ 0.0014 


it - 0.01256637 


- 0.0094 


- 0.0096 


+ 0.0002 


w - 0.01099557 


- 0.0951 


- 0.0960 


+ 0.0009 


7T - 0.009424777 


- 0.4270 


- 0.4300 


+ 0.0030 


tt - 0.0078539816 


- 0.4973 


- 0.4999 


+ 0.0026 


7T - 0.006283185 


- 0.1156 


- 0.1161 


+ 0.0005 


tt - 0.0047123889 


- 0.1260 


- 0.1267 


+ 0.0007 


tt - 0.00314159264 


- 1.6530 


- 1.6582 


+ 0.0052 


w - 0.0015707963 


- 6.4362 


- 6.4540 


+ 0.0178 


7T 


- 122.73 


- 58.54 




jt 4- 0.0015707963 


- 6.4362 


- 6.4184 


- 0.0178 


tt + 0.00314159264 


- 1.6530 


- 1.6481 


- 0.0049 


tt + 0.0047123889 


- 0.1260 


- 0.1255 


- 0.0005 


tt + 0.006283185 


- 0.1156 


- 0.1154 


- 0.0002 


tt + 0.0078539816 


- 0.4973 


- 0.4949 


- 0.0024 


* + 0.00942477 


- 0.4270 


- 0.4243 


- 0.0027 


tt + 0.01099557 


- 0.0951 


- 0.0944 


- 0.0007 


w + 0.01256637 


- 0.0094 


- 0.00964 


+ 0.0002 


t + 0.01413716 


- 0.1389 


- 0.1378 


- 0.0011 


tt + 0.0157079 


- 0.1662 


- 0.1646 


- 0.0016 
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Fig. 6 — Driven-wave loss versus A/3, near A/3\ B 
Aa/c | = 0,1,5. 



= 2tt (fundamental) for 



earlier perturbation theory 4 has shown that the fractional frequency 
interval between half-height points on the loss curve is 



A/ 
/ 



A/3 



*S 



A/3 



= 2 



Aa 

A/3 



(24) 



or 



A/3^ = 2 | Aa |. (25) 

Table VII shows a comparison between that limiting value and the 
true value for the numerical cases above, including | Aa/c | from one 
to 50. Even at I Aa/c I = 1 there is only a 30 per cent error. 
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Consider the limiting case Aa = 0. Then 

*i - [1 - TO 2 ]* 
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(26) 



and nulls in the loss curve occur when equation (19) representing 
E 2 is zero. This gives 



*. I...... -?V£)--i. 



(27) 



When we also have the perturbation condition, cz « 1, equation (27) 
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Fig. 7 — Driven-wave loss versus A/3, near A/3\ m = 2ir (fundamental) for 
Aa/c\ = 1,5,50. 
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becomes 



A3 



2tt 



* |null = 



(28) 



which agrees with previously known perturbation theory. 3 For this 

Table III — Exact and Approximate Calculations 
OF E x FOR I Aa/c I = 1 



Kt) 


20 log | Ei | 


Approximate 


Approximate 


Exact 


—exact 


it - 0.0157079 


- 0.2898 


- 0.2910db 


+ 0.0012db 


7T - 0.01413716 


- 0.3438 


- 0.3450 


+ 0.0012 


■k - 0.01256637 


- 0.4152 


- 0.4164 


+ 0.0012 


■k - 0.01099557 


- 0.5452 


- 0.5440 


+ 0.0012 


7T - 0.009424777 


- 0.7300 


- 0.7319 


+ 0.0019 


it - 0.0078539816 


- 0.9601 


- 0.9620 


+ 0.0019 


■k - 0.006283185 


- 1.3106 


- 1.3126 


+ 0.0020 


7T - 0.0047123889 


- 2.0216 


- 2.0243 


+ 0.0027 


7T - 0.00314159264 


- 3.3781 


- 3.3823 


+ 0.0042 


7T - 0.0015707963 


- 5.3475 


- 5.3523 


+ 0.0048 


it 


- 6.5701 


- 6.5705 


+ 0.0004 


ir + 0.0015707963 


- 5.3475 


- 5.3434 


+ 0.0041 


it + 0.00314159264 


- 3.3781 


- 3.3746 


- 0.0035 


7T + 0.0047123889 


- 2.0216 


- 2.0194 


- 0.0022 


ir + 0.006283185 


- 1.3106 


- 1.3092 


- 0.0014 


7T + 0.0078539816 


- 0.9601 


- 0.9588 


- 0.0013 


7T + 0.00942477 


- 0.7300 


- 0.7289 


- 0.0011 


w + 0.01099557 


- 0.5425 


- 0.5416 


- 0.0009 


tv + 0.01256637 


- 0.4152 


- 0.4146 


- 0.0006 


7t + 0.01413716 


- 0.3438 


- 0.3430 


- 0.0008 


7T + 0.0157079 


- 0.2898 


- 0.2892 


- 0.0006 



Table IV — Exact and Approximate Calculations 
of Ei for | Aa/c I = 5 



Afl/^ 


20 log 10 | Ex | 


Approximate 


*w 


Approximate 


Exact 


-exact 


w - 0.0314159262 
ir - 0.01884955 


- 0.2446db 

- 0.5337 


- 0.2452db 

- 0.5344 


+ 0.0006db 
+ 0.0007 


T 


- 1.6195 


- 1.6197 


+ 0.0002 


ir + 0.01884955 
ir + 0.0314159 


- 0.5337 

- 0.2446 


- 0.5344 

- 0.2445 


+ 0.0007 
+ 0.0001 
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Table V — Exact and Approximate Calculations 
of E x for | Aa/c | = 50 



Kt) 


20 logio I Ei | 


Approximate 


Approximate 


Exact 


-exact 


7t - 0.314159 
w - 0.1884955 

t + 0.1884955 
TV + 0.219911 


- 0.0234db 

- 0.0523 

- 0.1724 

- 0.0523 

- 0.0418 


- 0.0238db 

- 0.0525 

- 0.1727 

- 0.0525 

- 0.0421 


+ 0.0004db 
- 0.0002 
+ 0.0003 
+ 0.0002 
+ 0.0003 



case, from (26) 



E. ^ 1 - * I E 2 



= 1 



1 . 


sin 


x/M + i- 


'.' 




JW + i 



(29) 



and the loss has the form (just as in the case of sinusoidal coupling 3 

'sin au 



which in our terminology has half-peak loss (with A/3^ the variable) at 

1.8tt 



A0*j k-o = -^— 
For the numerical case of this paper, z = 1000 feet and 
A/3 #t Ua-o = 0.00565. 

Table VI — Phase Angle of E x for | Aa/c I = 



(30) 



Kt) 


Angle for E\ 


Approximate 


Approximate 


Exact 


-exact 


7t + 0.0015707963 
ir + 0.00314159264 
■k + 0.00471238891 

ir + 0.0157079 


- 67.753° 

- 47.111 

- 30.086 

- 9.085 


- 67.726° 

- 47.096 

- 30.092 

- 9.104 


- 0.027° 

- 0.015° 
+ 0.006° 

+ 0.019° 



1814 THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER li)(3S 



70 


\ 


















i 


















\ 




















W 
UJ 

a. 
o 

UJ 


" 


|Aa 
"1 c 


= 
















Z 


\ 


















u. 
O 30 


\ 






































UJ 

< 

I 

D_ 

20 










































/ \ 




















10 



/ |^&| = 1 X 


























z'' 




Aa 
C 


= 5 


— 





0.005 0.010 0.015 0.020 0.025 0.030 0.035 

A/9 # (2^\ IN RADIANS 



Fig. 8 — Phase of driven wave versus A/3.\ m /2 near Ai3X m = 2?r for | Aa/c | = 
0,1,5. 



In Table VII, the calculated value 0.0064 at | Aa/c | = 1 is reasonable, 
since D. T. Young's 4 work based on the perturbation theory indicates 
that the true loss peak is the convolution of the shape for | Aa/c \ » 1 
with the shape for | Aa/c | = 0. 

Consider now the peak loss at A/3 + = 0. When | Aa/c \ » 1, and 
e 4 "' « 1, it can be shown that equation (9) simplifies to 



Ei I = e" 



(31) 



Aa 
c 

A at 



»1 
«1 



tf* = o 
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where 



t 3 — -- 2 c 
x 



c 
Aa 



(32) 



In Table VIII we compare the loss computed from (32) with the 
actual loss, and see that even for | Aa/c | = 1 the error is only «30 
percent. A consideration of the terms of equation (9) indicates that these 
errors would be approximately constant with increasing z. 

A further simplification of the calculation of loss components now 
seems justified. For | Aa/c | ^ 1 it would appear that (31) and (32) can 
be used to calculate the peak loss due to a single square-wave coupling 
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Fig. 9 — Driven-wave loss versus A/3, near A/3Xm = 6jt (third harmonic) for 
Aa/c | = 0,1,5. 
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Fig. 10 — Driven-wave loss versus A/3, near A/3X m = 6ir (third harmonic) for 
| Aa/c | = 5,50. 

component. And the shape versus A/?* can be calculated from the work 
of D. T. Young, leading to an over-all attenuation defined by 



Ei I = e r " 



where 



r 4 = —zC 

7T 



C 

Aa 



i + if 1 



(33) 
(34) 



in which we require z and Aa/c such that e*"* « 1, | Aa/c | 2s 1. 
Table VII — Limiting and True Values 



(1) 

1 A " 

1 c 1 


(2) 

Limiting value of A/S.i/j 
from equation (25) 


(3) 

Calculated A0.i/i 
from Figs. 6 and 7 


Ratio of 

Col. (3) 

Col. (2) 


1 

5 
50 


0.00493 

0.0246 

0.246 


0.0064 
0.0265 
0.250 


1.3 

1.077 

1.016 
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Table VIII — Loss Comparison 



A)3. = 




Aa 
C 




Loss (dB) from 
equation (32) 


Actual loss (dB) 
equation (9) 


1 

5 
50 


- 8.68 

- 1.74 

- 0.1737 


- 6.57 

- 1.62 

- 0.1727 
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appendix a 



Coupling Functions 

We use the following equations to represent two coupled waves of 
amplitude E-l and E 2 : 



Et(z) = -y 1 E 1 + c 21 (z)E 2 



j z E 2 {z) = e 12 {z)E 1 -y 2 E 2 

where y„ = a n + j(3 n = complex propagation constant. Let 

c»i(«) = ic P e' kcZ . 
Energy conservation leads to: 

c.o(z) = jc v e~ ik ". 
Also let 

E x = e-^'V, 
E 2 = e-"'V 2 . 
Then (35) and (36) become: 

dz JCp - e 
dz JCpVie 



(35) 
(36) 

(37) 
(38) 

(39) 

(40) 

(41) 

(42) 
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We now see that the solution for k c = holds for k nonzero provided 
that we make the change 

AB = ft - ft + k c 

and put AB in place of Aft = ft — /? 2 in the solution for k B = 0. 

Complete transfer of power between waves can occur when AB = 
or 



K = (ft - ft). 



(43) 



APPENDIX B 



Square-Wave Coupling 

We start with equations (1) through (4) representing uniformly 
coupled waves with k e = (see Fig. 1). Equations (9) through (14) with 
c* = c and Aft. = A/3 give the amplitude of the output waves at z for 
the conditions E x = 1.0 and E 2 = at z = 0. We apply these equations 
to the coupling distribution of Fig. 2 and find that the output amplitude 
for the undriven wave at z = X m /2 is: 



E 2 | 2=Xm /2 = exp 



-iifix + ft) f 



} 



+isin 



cX, 



•^ + 89 



v» + (M 



(44) 



In the above, the simplification Aa = has been assumed. 

Using the output waves at z = X,„/2 as input conditions to the 
following coupling region (Fig. 2), the amplitude of the undriven 
wave at z = X m is found to be: 



E, 



= (-D 



i§ 



[& - • 



k2c I 



The ratio of the undriven wave amplitude at z = A™ to the undriven 
wave amplitude at z = Xm/2 is [from (44) and (45) ] 
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E 2 L . *\2c 

Hi 9 X_/5 



Vi + « 



"><-Nl + (ff,_ 



exp 



; (ft + ft) 

■3 a K 



(46) 



By imposing the condition 



we notice that the magnitude of (46) approaches two for A/3 }>> 2c, that 
is, for very small coupling. Fig. 4 sketches the undriven wave amplitude 
as a function of z during the first two coupling intervals. We will now 
express the approximate or average coupling between the waves by 
using the linear approximation E 2 (z) = 2cz/t. We notice that uniform 
coupling without phase reversal would have resulted in the relation 
E 2 (z) = cz. We therefore arrive at the transformation 

c (constant coupling) becomes 

2c/ ir (for periodically reversed coupling). 

(48) 

The associated transformation of the condition for maxmium energy 
transfer, from (47) , is 



We might notice here that the departure of the actual amplitude in 
Fig. 4 from the straight-line approximation, shown by e in Fig. 4, has a 
maximum value which can be shown to be 0.21(cX m )/7r. Thus the devia- 
tion between our straight-line approximation and the actual amplitude 
becomes smaller for diminished values of coupling per unit X m . 

We now specify c^ and A/3^ in equations (9) through (14) to represent 
the waves with periodically reversed coupling using the average coupling 
approach. Since the in-phase build-up of power in the undriven wave is 
a maximum for A/3 specified by equation (49), we define a new dif- 
ferential phase parameter to give a departure from this condition: 



V.-6fi-*y]l -($*)'■ 
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By comparison with equation (7), which represents an exact solution 
for complex coupling, this definition of A/3* seems reasonable. Numerical 
checks reported in the text verify this presumption. The value of c* is 
given by (48). 

It has already been noted, near equation (46) that two successive 
X m /2 coupling intervals give twice the undriven line amplitude com- 
pared with the first X m /2 interval, which is just what occurs for small 
constant coupling (no reversal) and A/3 = 0. We also find that A/3* and 
c* correctly give the wave amplitudes for periodically reversed coupling 
at 2 equal to integral multiples of X m /2, even when A/3 is not large com- 
pared with 2c. For example, letting E 2 at z = X m /2 be 0.707 requires 
(A/3/2c) = 1 (see equation 44) ; we maintain equation (47) and (48) as 
before. Then, from (45), | E 2 | becomes unity at z = X m which is also 
predicted by A/3* and c* in (9) through (14) and which is analogous to 
the behavior of two conventional 3 dB directional couplers in cascade. 

The above discussion represents the changes in wave propagation 
introduced by coupling for A/3 in the vicinity of the value given by 
equation (49). 

The perturbation solution 3 for sinusoidal periodic coupling is known 
to yield coupling absorption peaks when 

A/3X m = 2irp ^ 

V - 1, 3, 5 ••• . 

Similarly, there are other regions of strong interaction for the square- 
wave coupling of Fig. 2. For example, consider Fig. 5, which represents 
the situation when 



kc\ m _ Sir 
2 ~ 2 ' 

There is another region of A/3 defined by 



(52) 



A/3L„ imuroconver3ioa = 3^Ji-^) (M) 



K \ x V37T 

where there is a local maximum of conversion. As diagramed in Fig. 5, 
the average conversion coefficient is 2c/37r. Thus the appropriate values 
of A/3* and c* for equations (9) through (14) are 

:l77 (54) 



'* 



A/3 = A/3-3- ^Jl (<*■ 



^V 1 " V3x 



A/3* - 


2c 

c = — 

* pir 


xpir 1 
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More generally there are absorption peaks at 

(55) 

(56) 
where p = 1, 3, 5, . . . 

APPENDIX C 

Sine Wave Coupling 

We start with equations (1) and (2) with the coupling defined 

c 2 i = c 12 = jc sin ( — j • (57) 

Using perturbation theory 3 , and letting A« = 0, 

Bl = c f sin (~V '*" ds (58) 



which yields 



E, = _, ,_" ^ <;- - «"' A ' 



m - 



A/3' 



ft -© + "* &)]}■ 






Evaluating at A/? = 27r/A„, and 2 = nA„,/2 yields 

ftu^-.-A. - -4%" (60) 

with n = 1,2,3.... 

It can be verified that A/3 = 2w/\„, yields the maximum value of 

E-2 at z — nkm/2. 

Thus the equivalent uniform coupling value for sinusoidal coupling is 

c, = c/2 (61) 

which appears in Table I. 
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